# -*- mode: snippet -*-
# name: lme from nlme Linear Mixed-Effects Models
# key: nlme.lme
# contributor : Shuguang Sun
# --
res.$1 <- nlme::lme(fixed  = ,
              data   = $1,
              random = ~ | id,
              method = "REML"
              )

summary(res.$1)
intervals(res.$1)
getVarCov(res.$1)
